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Abstract 

We study the Hawking radiation in the framework of open quantum systems by examining the 
time evolution of a detector (modelled by a two-level atom) interacting with vacuum massless 
scalar helds. The dynamics of the detector is governed by a master equation obtained by tracing 
over the field degrees of freedom from the complete system. The nonunitary effects are studied 
by analyzing the time behavior of a particular observable of the detector, i.e., its admissible state, 
in the Unruh, Hartle-Hawking, as well as Boulware vacua outside a Schwarzschild black hole. We 
find that the detector in both the Unruh and Hartle-Hawking vacua would spontaneously excite 
with a nonvanishing probability the same as what one would obtain if there is thermal radiation 
at the Hawking temperature from the black hole, thus reproducing the basic results concerning the 
Hawking effect in the framework of open quantum systems. 

PACS numbers: 04.70.Dy, 03.65.Yz, 04.62. +v 
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I. INTRODUCTION 



Black holes are intriguing objects and are worth studying in all possible varieties. The idea 
of black holes has been proven to be highly fruitful. In particular, Hawking's discovery that 
black holes are not, after all, completely black, but quantum mechanically, emit radiation 
with a thermal spectrum [1], has provided us with the understanding that black holes may 
play the role of "Rosetta stone" to relate gravity, quantum theory and thermodynamics. 
Therefore, Hawking radiation, as one of the most striking effects that arise from the combi- 
nation of quantum theory and general relativity, has attracted widespread interest in physics 
community and it has been extensively examined from different prospectives, yielding dif- 
ferent derivations of it. These derivations include (but not limited to) Hawking's original 
one which calculates the Bogoliubov coefficients between the quantum scalar field modes of 
the in vacuum states and those of the out vacuum [l, 2], an Euclidean quantum gravity 
derivation [3] which has been interpreted as a calculation of tunnelling through classically 
forbidden trajectory [1], an approach based upon string theory [5, 6], an interesting proposal 
which ties its existence to the cancellation of gravitational anomalies at the horizon [7] , and 
a recent study which reveals an interesting relationship between the existence of Hawking 
radiation and the spontaneous excitation of atoms using the DDC formalism [8] that sepa- 
rates the contributions of vacuum fluctuations and radiation reaction to the rate of change 
of the mean atomic energy [9] . 

In the current paper, we shall try to understand the Hawking radiation by examining the 
time behavior of a static detector (modelled by a two-level atom) outside a Schwarzschild 
black hole immersed in vacuum massless scalar fields using the well-known techniques devel- 
oped in the study of open quantum systems. As for any open system, the full dynamics of 
the detector can be obtained from the complete time evolution describing the total system 
(detector plus external fields) by integrating over the field degrees of freedom, which are in 
fact not observed. It is worth noting here that an examination of a similar issue, i.e., the 
Unruh effect associated with uniformly accelerated atoms in the paradigm of open quantum 
system has been already been carried out [10]. 
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The paper is organized as follows. In next Section, we shall review the basic formalism, 
the derivation of the master equation describing system of the detector plus external vacuum 
scalar fields in weak coupling limit and the reduced dynamics it generates for the finite time 
evolution of the detector. In Section III, we apply the method and results of the preceding 
Section to discuss the probability of spontaneous transition of detector from the ground 
state to the excited states outside a Schwarzschild black hole. Finally, we conclude with 
some discussions in Section IV 

II. THE MASTER EQUATION 

We shall consider the evolution in the proper time of a static detector (two-level atom) 
interacting with vacuum massless scalar fields outside a Schwarzschild black hole and assume 
the combined system (detector + external vacuum fields) to be initially prepared in a factor- 
ized state, with the detector held static in the exterior region of the black hole and the fields 
in their vacuum states. Our derivation of the master equation in this Section follows closely 
to that in Ref. [10]. The atom is assumed to be fully described in terms of a two-dimensional 
Hilbert space, so that its states can be represented by a 2 x 2 density matrix p, which is 
Hermitian p + = p, and normalized Tr(p) = 1 with det(p) > . In order to achieve a rigorous, 
mathematically sound, derivation of the reduced dynamics of the detector, we shall assume 
that the interaction between the detector and the scalar fields are weak so that the finite- 
time evolution describing the dynamics of the detector takes the form of a one-parameter 
semigroup of completely positive maps [11, 12]. Without loss of generality, we take the total 
Hamiltonian for the complete system to have the form 

H = H S + H (j> + \H' . (1) 

Here H s is the Hamiltonian of the atom, which in the most generic case takes the form 

3 

1=1 

where <7j (i = 1,2,3) are the Pauli matrices, ojq the energy level spacing and n = (m, n^, n 3 ) 
a unit vector. In the present paper, we will take, for simplicity, n to be along the third axis 
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such that H s simplifies to 



H s = f . (3) 



Hfj, is the standard Hamiltonian of massless, free scalar fields, details of which need not be 
specified here and H' is the interaction Hamiltonian of the atom with the external scalar 
fields and is assumed to be given by 

H' = a^{x) . (4) 

It should be pointed out that the coupling constant A in (1) is small, and this is required by 
our assumption that the interaction of the atom with the scalar fields is weak. 

Initially, the complete system is described by the total density p tot = p(0) <g> |0)(0| , where 
p(0) is the initial reduced density matrix of the atom, and |0) is vacuum state of the field 
In the frame of the atom, the evolution in the proper time r of the total density 
matrix, p tot , of the complete system satisfies 

= -iL H [ P Ur)\ , (5) 

where the symbol Lh represents the Liouville operator associated with H 

L H [S] = [H,S}. (6) 

The dynamics of the atom can be obtained by tracing over the field degrees of freedom, i. 

e., by applying the trace projection to the total density matrix p(r) = Tr$[p tot (r)] . 

In the limit of weak-coupling which we assume in the present paper, the reduced density 

is found to obey an equation in the Kossakowski-Lindblad form [13, 14] 

d P (r) 
dr 

where 



: [H eS: p(r)] + C[p(r)] , (7) 



£>[P] = ^ a ij [ 2 a jP °i - °i a 3 P~ P a i a j] ■ ( 8 ) 

The matrix and the effective Hamiltonian H c r are determined by the Fourier and Hilbert 
transforms of the field vacuum correlation functions 

G + (x-y) = (0mx)$(y)\0), (9) 



which are defined as 



0(A) = 
/C(A) 



dre lXT G + (x(t)) 

-oo 

7ri 



du 



— A 



Then the coefficients of the Kossakowski matrix be written explicitly as 



(10) 



'111 



(12) 



with 



A = l[g(u ) + g(-w )] , B = l[Q{u> Q ) - g(-co )} , C = 0(0) - A . (13) 



Meanwhile, the effective Hamiltonian, H e s, contains a correction term, the so-called Lamb 
shift, and one can show that it can be obtained by replacing in H s with a renormalized 
energy level spacing Q as follows [10] 



H cS = ^0-3 = {^0 + i[/C(-w ) - /C(cjo)]}o-3 



(14) 



where a suitable substraction is assumed in the definition of K{—ojq) — K{uq) to remove the 
logarithmic divergence which would otherwise be present. 

To facilitate the discussion of the properties of solutions for (7) and (8), let us express the 
density matrix p in terms of the Pauli matrices, 

P(r) = + Y^Pi(T)a^ . (15) 

Substituting Eq. (15) into Eq. (8), it is easy to show that the Bloch vector \p(r)) of compo- 
nents {pxi^r) , p 2 (r) , p 3 (r)} satisfies 





Or 



\p(r)) = -2H\p(r)) + \ V ) , 



(16) 



where \rj) denotes a constant vector {0, 0, —45}. The exact form of the matrix Ti reads 



( 2A + C Sl/2 \ 



n 



-Q/2 2A + C 
2A 



(17) 



Eq. (16) can be solved formally, with the solution being written as 

|p(r)) = e- 2 ^|p(0)) + (l-e" 2 ^)| Poo ), 



with 



IP-) = \n- l \q) = ~ 





w 



(19) 



Here the matrix (operator) e~ 2Hr is defined by series expansion as usual. However, since H, 
obeys a cubic eigenvalue equation, so powers of TC higher than 2 can always be written in 
terms of combinations of 7i 2 , 7i and /, the 3x3 unit matrix. Actually, three eigenvalues of 
7i can be found explicitly: \\ = 2 A, A± = (2 A + C) ± iQ/2. One can then show that 



-2Hi 



4 



tt 2 + AC 2 



e-^A, + 2e~ 2( - 2A+c > 



A 2 cos(fir) + A ; 



sin(Or) 



(20) 



where 



n 2 

Ax = [{2A + C) 2 + — ] I - 2{2A + C)H 



A, 



-2A(A + C)I + (2A + c)n - -n 2 



n 2 n 2 

A 3 = 2A[— - C(2A + C)]I+ [C{AA + C) - — ]H - CH 2 



(21) 



Eq. (20) reveals the effects of decoherence and dissipation on the atom characterized by the 
exponentially decaying factors involving the real parts of the eigenvalues of 7i and oscillating 
terms associated with the imaginary part. These nonunitary effects can be analyzed by 
examining the evolution behavior in time of suitable atom observables. For any observable 
of the atom represented by a Hermitian operator O, the time behavior of its mean value is 
determined by 

(0(t)) = Tr[Op(r)} (22) 

Let the observable O be an admissible atom state pf, then Eq. (22) yields the probability, 
Vi-tffr), that the atom evolves to the expected state represented by density matrix Pf(r) 
from an initial one p(0) = 

V^ f (r) = Tr[ Pf p(r)} . (23) 
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If initially the atom is in the ground state, so that its Bloch vector |p(0)) is {0, 0, —1}, and 
the final state pf is the excited state given by the Bloch vector \pf) = {0, 0, 1}, then one has, 
by making use of Eqs. (18-23), 



1 f. ^ Ar \ (\ B 



p ww=2i'- e jrij (24) 

The probability per unit time of the transition from the ground state to the excited state, 
in the limit of infinitely slow switching on and off the atom-field interaction, i.e., the spon- 
taneous excitation rate can be calculated by taking the time derivative of Vi^>f(r) at r = 0, 



= 2A-2B = 2g(-uj ) . (25) 

T=0 



III. PROBABILITY OF SPONTANEOUS TRANSITION OF THE DETECTOR 
OUTSIDE A SCHWARZSCHILD BLACK HOLE 

Now, we apply the open quantum system formalism developed in the preceding Section to 
address the issue of finite time evolution of a static detector (two-level atom) interacting with 
vacuum scalar fields outside a spherically symmetric black hole and calculate the probability 
of spontaneous transition and the spontaneous excitation rate of the detector from the ground 
state to the excited state. The metric of a spherically symmetric black hole can be written, 
in the Schwarzschild coordinates, as 

ds 2 =(l- —) dt 2 - - tl/ " r 2 (d9 2 + sin2 6d <t> 2 ) ■ ( 26 ) 



r J 1 - 2M/r 

However, a delicate issue that we have to address first before moving on is how the vacuum 
state of the quantum fields is determined. Normally, a vacuum state is associated with 
non-occupation of positive frequency modes. However, the positive frequency of field modes 
are defined with respect to the time coordinate. Therefore, to define positive frequency, one 
has to first specify a definition of time. In a spherically symmetric black hole background, 
one definition is the Schwarzschild time, t. The vacuum state, defined by requiring normal 
modes to be positive frequency with respect to the Killing vector d/dt with respect to which 
the exterior region is static, is called the Boulware vacuum [15]. Other possibilities that have 



been proposed are the Unruh vacuum [16] and the Hartle-Hawking vacuum [17]. The Unruh 
vacuum is defined by taking modes that are incoming from to be positive frequency 
with respect to d/dt, while those that emanate from the past horizon are taken to be positive 
frequency with respect to the Kruskal coordinate u, the canonical affine parameter on the 
past horizon. The Hartle-Hawking vacuum, on the other hand, is defined by taking the 
incoming modes to be positive frequency with respect to v, the canonical affine parameter 
on the future horizon, and outgoing modes to be positive frequency with respect to u. 
Therefore, in what follows, we shall calculate the probability of spontaneous transition of a 
static detector from the ground state to its excited state and spontaneous excitation rate in 
these three vacuum states. 



A. The Unruh vacuum 



We start with the Unruh vacuum. The Wightman function for massless scalar fields in 
the Unruh vacuum can be computed as follows [18, 19, 20] 



Gu + (x,x') = y 



Yr 



lm\ 



I Ri(u>,r) 



(27) 



where k = 1/(4M) is the surface gravity of the black hole. The Fourier transform with 
respect to the proper time, which is related with the coordinate time through, 



dr = Jg^dt 



2M , 
1 dt , 



(28) 



reads 



e lXT Gu + {x,x')dT 



1 oo 
1=0 



8(\y/g^)(l + 2l)\R l (\ y /g^,r 



[l + 2l)\R l (X y /g^,r] 

I _ g-27rA v /goo/K 



|2l 



(29) 



where we have used the relation 



y \ Y lm (e,, 



m=—l 



21 + 1 

47T 



(30) 



The summation in Eq. (29) is not easy to evaluate into a closed form in general. However, 
its behavior both close to the event horizon and in the spatial asymptotic region, which, 
fortunately, are regions we are most interested in, can be found using the following properties 
of the radial functions 

4u 2 



^(2/+i)i^(^or~ < 



I — 2M. ' 



1=0 



J2(2l + l)\%(cu,r) 



1 oo 

-^(2/ + l)|^M| 2 

1=0 
1 oo 

_£ (a+ l)|B, H 



z=o 



and 



E 

1=0 



Auj 2 



B,(u) 



r ->2M 



r — > oo 



r — »• oo 



(31) 



10 



4M5 



(32) 



(33) 



uj^O 



With the above expression, it is easy to show that 



iirujQ 



4wn 



_ 9_M 

\ — e -2iTLu /K r 4M 2 



£)(l + 20|£i(aWffM) 



1 + COth (^ /K r ) - + ^g^^ |2 + WM 
^— ' 27T 



167ro;or 2 



r -> 2M 



r — > oo 



(34) 



where K r = k/v/1 — 2M/r . This leads to the following behaviors of the coefficients of the 



Kossakowski matrix 



r -> 2M : < 



oo 



B 



v 



oo 



64M 2 7ru/ 



C[7 = MO) - Ay . 



/ — ^ i2 , w o 



(35) 



r — > oo : < 





^0 

~ T~ + 


F (+ Vo) 


47T 




Bu 


^0 . 




47T 




C v 


= Mo) 


-At/ , 



(36) 



where the auxiliary function is defined as 
-1 + coth(^a) 



u ) 



327ru;r 



E( 1 + 2/ ) 



1=0 



(37) 



Substitute the above coefficients into Eq. (24) and Eq. (25), then the probability of spon- 
taneous transition to the excited state of the ground state detector and the spontaneous 
excitation rate can be obtained respectively 



J) 



-^fQ-coth nrwo//Cr) 



1 



+ 



#00 (e 



1) 



and 



1=0 



I _|_ e 2iru) /K r 2(1 + g27rwo/Kr) 

2/)| J B z (-^ 0v ^)| 2 



2M 



OG 



u 



2M 



oo 



(38) 



yj- I K r 

/(-wb,r)w 



where 



f(u ,r) = 



l + 2l)\B l (u 0y /tod | 2 



z=o 



4u;n 2 r 2 



(39) 



(40) 



The above results reveal that, close to the horizon, the ground state detector in the vacuum 
would spontaneously excite with an excitation rate same as what one would expect if there 
is a flux of thermal radiation at the temperature T = k t /(2tt). However, as the atom is 
placed away from the horizon, this thermal flux is backscattered by the spacetime curvature, 
resulting in depletion of part the flux, the effect of which is described by the greybody factor 
f(—u>o, r). The depletion is dependent on the atom's radial distance from the black hole and 
becomes greater as the detector is placed farther away. The effective temperature T, which 
can be understood as a result of combined effects of thermal radiation from the black hole 
and the Unruh effect which results from the acceleration with respect to the local free-falling 
inertial frame needed to hold the detector static at a finite radial distance, approaches 
the Hawking temperature T H = n/(2ir) in the spatial asymptotic region, suggesting the 
temperature of the thermal radiation emanating from the horizon is just the Hawking value, 
since the acceleration needed to hold the detector static vanishes in the spatial infinity This 
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understanding we gain for the Unruh vacuum in the paradigm of open quantum systems is 
consistent with what we have in other different contexts [20, 21, 22]. 



B. The Hartle-Hawking vacuum 



Now, let us turn to the Hartle-Hawking vacuum case. Then the Wightman function for 
the scalar field becomes [18, 19, 20] 



G H + {x,x) = ^ 



ml 



+ 



°° I Yi m ( 



A'KUJ 



The Fourier transform is given by 

POO 

g H (X) = / e iXT G H +(x,x')dr 



Ri(u,r) 



g— iujAt ^ 

— 7-I RAuj, r) 



(41) 



E 



;i+2Q 

8ttA 



Ri{^y/<toa,r)\ 2 , l^(-^v / ^00' r 



12 n 



(42) 



which can be evaluated approximately at the event horizon and in the spatial asymptotic 
region to get 



' 1 + coth (^a) 
16^0 

1 + coth (if) 



w + 



;i + 2/)|A(-^ 0v ^ 



Z=0 



AM 2 



, r -> 2M , 



E 



:i + 20|5,(woV5ro) 



+ 4cj #00 



(43) 



r — ► oo 



The corresponding coefficients for then read 

cj coth(^a) -l + coth(^a) 



A 



ii 



Air 



^TTUlQ I K r I g / 



128M 2 7rcu 



E( 1 + 2/ ) 



«=o 



r^2M: < 



cj n -l + coth(^a)^ 



4tt 128M 2 7tcj 

+ e W«r| Sl (_ a;()V ^ 

Cif = fe(0) - Atf , 



i=0 



(44) 
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^~_ coth (__) + _ _ 



oo : < 



Air k r z 

ujq_ F^jup) (45) 



B H pa — , 

47r r z 



k cw = mo) - a h • 

Therefore, the probability of spontaneous transition and the spontaneous excitation rate 
from ground state to exited state are respectively 



g-^ 2 - coth(7rajo/re r ) 

1 _|_ ^TTUlo/Kr 



2M 



and 



1 _|_ e 2-KU(_)/K 

u f{u> , 2M) uj 



OO 



2M 



w | /(-w ,r)w 

"T T — ^ OC 

Here, one can see that the spontaneous excitation rate consists of two pieces. One is the 
thermal radiation outgoing from the horizon (the first term in when r —>■ 2M), which 
is identical to what we have obtained in the Unruh case, and the other can be viewed as the 
thermal radiation incoming from infinity (the first term in when r — > oo). Both the 

outgoing and incoming thermal radiation are backscattered by the curvature on their way 
to a finite radial distance from the horizon, and this backscattering is represented by the 
greybody factors in the second term in for both cases. Therefore, one sees that the 
Hartle-Hawking vacuum actually corresponds to a black hole immersed in a bath of both 
incoming and outgoing thermal radiation. Again, this is consistent with our understanding 
of the Hartle-Hawking vacuum gained in other different contexts. 



C. The Boulware vacuum 

Finally, let us briefly discuss the Boulware vacuum case. Now, the Wightman function 
can be written as [18, 19, 20] 

^-—| Y lm (9,</>) I 2 [I Ri(u,r) | 2 + | R^u.r) \ 2 }du , (48) 

im 
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The corresponding Fourier transform reads 

/oo 
e iXT G B + [x(r)}dr 
-oo 

= E 



mi 



du 



i(A-u/ v /l-2M/r)T| y 



Ri(uj,r) 



ml J0 



u/^l-2M/r)\Y l 



lm \ 



Ri(uj,r) 



\Ri(uJ,r)\ 2 ]dT 
Ri(u,r)\ 2 } .(49) 



The further calculation should proceed in three separate cases 



Qb{\) = { 



ml 
ml 



9oo, r 



Ri(XV9Vo,r)\ 2 ] 



Ri(*V%o~, r)\ 2 + \ %(Xy/g^, r) | 2 ] 



A^O 



0. 



A > , 

A = , (50) 
A < . 



This indicates that for a positive Uq, Qb{— ^o) vanishes. Thus a substitution of Eq. (50) 
into Eq. (13) yields that As = Bb , which in turn gives rise to Vi^f(r) = and = 
2Ab — 2Bb = . Therefore, no spontaneous excitation would ever occur in the Boulware 
vacuum. As a result, the Boulware vacuum corresponds to our familiar notion of a vacuum 
state. 



IV. CONCLUSION 

In summary, we have examined the Hawking radiation from the point of view of open 
quantum systems by looking at the time evolution of a detector (modelled by a two-level 
atom) interacting with vacuum massless scalar fields. The time evolution of the detector is 
governed by a master equation obtained by tracing over the field degrees of freedom from the 
complete system. The nonunitary effects have been studied by analyzing the time behavior 
of the detector's admissible state in the Unruh, Hartle-Hawking, as well as Boulware vacua 
outside a Schwarzschild black hole. It is found that the detector in the Unruh and Hartle- 
Hawking vacua would spontaneously excite with a nonvanishing probability the same as 
what one would obtain if there is thermal radiation at the Hawking temperature from the 
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black hole, reproducing the basic results concerning the Hawking effect in the framework of 
open quantum systems. However, as we have pointed out in Section II, the possibility of 
a nonvanishing spontaneous excitation rate in fact involves the phenomena of decoherence 
and disspation, so the open system approach seems more comprehensive physically than 
traditional treatments. Our study suggests that the general techniques and results in the 
theory of open quantum systems not only is applicable to the study of the Hawking effect and 
possibly even to that of other phenomena in curved space-times, such as particle creation, 
but also may shed new light on the physical understanding of them. We hope to turn to 
these issues in the future. 
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